SU(3) Quantum Interferometry with single-photon input pulses 
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We develop a framework for solving the action of a three-channel passive optical interferometer 
on single-photon pulse inputs to each channel using SU(3) group-theoretic methods, which can be 
readily generalized to higher-order photon-coincidence experiments. We show that features of the 
coincidence plots vs relative time delays of photons yield information about permanents, immanants, 
and determinants of the interferometer SU(3) matrix. 

PACS numbers: 



An optical interferometer is a coherent scatterer of 
light comprising passive optical elements such as phase 
shifters, mirrors, and beam splitters. The interferometer 
operates in the quantum regime by injecting nonclassi- 
cal light into the input ports and yields nonclassical light 
from the output ports. The most famous quantum op- 
tical effect is the two-photon Hong-Ou-Mandel (HOM) 
dip [l|, with the "dip" corresponding to extinction of 
photon-coincidences at the two output ports of a bal- 
anced beam splitter given two identical input photons 
at each input ports. Instead the two photons exit in a 
superposition state of leaving together from each of the 
two output ports. Important applications include charac- 
terizing distinguishability between pairs of independent 
photons measuring coherence Q and purity Q of 
single photons, producing two-photon entanglement [f|, 
performing dense coding Q and singlc-qubit quantum 
fingerprinting Q , and creating non-deterministic nonlin- 
ear gates in optical quantum computing 0] . 

Generalizing to higher-order photon-coincidence dips 
will be valuable in many ways including determining 
distinguishability between multiple photons simultane- 
ously, computing the permanent of special unitary (SU) 
matrices through photon coincidence dips Q, and effi- 
ciently sampling permanents of sub-matrices of SU ma- 
trices, which would be hard classically (with complex- 
ity BPP# p )[l3l. Experimentally eight individual pho- 
tons have been manipulated with exquisite control 
and three photons in three coupled interferometers have 
shown nonclassical interference 12| so higher-order pho- 
ton coincidence dips are feasible well beyond the two- 
photon two-channel case. Although coherent scattering 
of Fock states has been analyzed theoretically for vari- 
ous instances of interferometers [l3l - fl5j , a full multimode 
analysis of multi-channel interferometry incorporating re- 
alistic source and detector spectral responses and allow- 
ing for arbitrary configurations has not yet been fully 
studied despite its necessity for quantitative studies of 
these nonclassical interferometric systems. 

Here we use group theory to solve the general prob- 



lem of photon coincidence probabilities at the output of 
any passive three-channel interferometer (whose action is 
represented by the SU(3) matrix) acting on single-photon 
pulse inputs, and we show specifically how the three- 
photon coincidence dip can be generated. We exploit 
the fact that an input of at most one photon per mode 
yields an output state that is a superposition of states of 
irreducible representations (irreps) of SU(3), and we find 
that, conditioned on one photon per mode in the input 
and output with frequencies unpermuted, the amplitude 
of the transformation equals a sum of the permanent, 
immanants and determinant of the SU(3) matrix. Then 
extinguishing this weighted sum is a necessary and suffi- 
cient condition for a three-photon dip. 

To illustrate we construct one instance of a SU(3) inter- 
ferometer coincidence-rate dip by extinguishing the coin- 
cidence corresponding to one photon exiting each of the 
three output ports by judiciously choosing the eight gen- 
eralized Eider angles. Our results suggest a rich group- 
theoretic structure to photon interferometry and a direct 
relationship to permutation symmetries of matrices. Fur- 
thermore our methods create a pathway to solving for 
arbitrary n-mode interferometers rather than restricting 
to three modes as reported here. 

A single photon pulse with light-source spectral func- 
tion (j> (uj) in one spatial mode, or channel, is 



dc^ s H|iM), / da^ s (w)| 
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for | l(u>)) = a'(uj) |0) with a(uj) the creation operator at 
frequency ui satisfying [a,k((^i), aj (uij)] = 5kiS(uji — 
k and I labels for distinct channels, and 1 the identity 
operator. The parenthetical (rounded) bra-ket notation 
distinguishes frequency-explicit states from other states. 
For convenience, we choose a Gaussian spectral functions 
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cxp[-( W - Wo ) 2 /(4<7 2 )], (2) 



but our approach accommodates any physical spectral 
function. 
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The detector has a response function (jp(uj), possi- 
bly different for each output port. For convenience we 
take D (w) to be Gaussian, similar to ©, but perhaps 
with different mean and variance; the variance is in- 
versely related to the detector integration time. For 

f[J d^0 D K)at(^-)j |0) (3) 

an ?i-photon state corresponding to a superposition 
of single-photon states in different infinitesimal fre- 
quency modes, the ideal detector executes projective 
measurements of the type II„ := |n)(n|, J2^=o^ n 
i . Detector inefficiency arises due to a spectral mis- 
match between source and detector spectral functions: 
the probability of detecting a single source photon is 

Jdio[4> s {io)]*4> D {io) 2 < l. 

We assume identical photons from the source, with dis- 
tinguishability introduced by time delays of the pulses. 
A tunable time delay of r on the single-photon pulse 
is expressed in the Fourier domain by a phase shift 
|1) h-> J dw</>(w) exp[— ilut]\ 1(w)), and the inverse Fourier 
transform of </>(w)exp[— iur] is 4>{t — r). 

Now that we have established the single-mode source 
and detector formalism, we consider two photons imping- 
ing on the two input ports of the passive two-channel 
SU(2) interferometer resulting in a photon-coincidence 
dip. This transformation can be expressed as 
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with f2 comprising the three Euler angles a, /?, 7, and 
the optical-element transformations are assumed to be 
independent of frequency. 

The input state |11) S has one photon in each port, i.e., 
a tensor product of two single-photon source states ([T]) 
with time delay r between modes 1 and 2. As we 
need to convert the direct product of irreps for single- 
photon interferometer inputs to a direct-sum decomposi- 
tion, Young tableaux are valuable, especially for SU(n) 
interferometry with n > 2 as these cases can be compli- 
cated to construct. Young tableaux are diagrams com- 
prising Ds arranged in rows and columns that represent 
various partitions of n. 

The representation of any SU(2) matrix on two pho- 
tons with frequencies w\ and W2 decomposes as 



□ ® □ -i- m © □ 
(1,0)® (1,0) ->(2, 0)0 (0,0) 



(5) 



with the second row giving the standard labels of the 
irreps in terms of highest weight. The decomposition is 
also written as 3 © 1 with I = 0, 1 irreps (singlet and 
triplet, respectively) and m = —I, ...,£. 



The two-photon state is a combination of \lm) states: 
I l(wi)l(w2)) = ^ (|00) + |10». The single-photon state 

a\(u>i) |0) transforms under SU(2) as \t = 1 1 
a\{tjJi) |0) transforms as | — \ — |}.: 
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R{$l)a\(u x ) |0) = \a\{u x )D\ ^Q) + 8^{u 1 )D» 1 .(Q)} |0) 

L 2 ■■ 2 2 ■ 2 J 

R(n)at(uj2) |0) 

= a\(uJ2)Dj ^.(ty + aU^DK _ i (0)l |0) (6) 

2 ' 2 2 ' 2 J 

with D l m , m (Q) := (£m'\R(n)\£m) an element of the 
(2£ + l)-dimensional Wigner D-matrix. The two-photon 
transformation is thus 

R(Q)\ 1( W1 )1(W2)) = ^ (JJ(fi) [00) + R(Q) |10>) . (7) 

Extinguishing coincidences arising from the triplet 1 10) 
contribution, which is symmetric under an tui <-> L02 
exchange, requires that (10| R(Q) |11) = -^Dq (£1) = 
cos f3 = 0. Therefore, (3 = ^ir and is independent of a, 7. 
These parameters correspond to a balanced beam split- 
ter with a relative phase shift between modes 1 and 2. 
For a = 7r/4 = —7, the beam splitter is balanced and 
symmetric: B := i?(f2) for this important case. 

Typically detectors at different output ports are dis- 
similar. If the source pulse has carrier frequency ojo and 
bandwidth <jq, and detectors have carrier frequencies uii 
and bandwidths <Ji for each output mode i, we introduce 
weighted variances a 2 i := <Jq + af, 
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a a +(T t \<7 a, 

and Gaussian spectral mismatch function 

/ (ujq - U!i) 2 
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The two-mode output coincidence rate is 

P n = s (ll|S t ni (gillxBlllf 
=AiA 2 e ^ r2 /^i e ^i _ 



(8) 

(9) 
(10) 



Pn = as expected for indistinguishable photons and 
identical detectors. 

Now we proceed to three-channel SU(3) interferometry 
transformation, with the eight-parameter generalized Eu- 
ler angle f2 (assumed to be frequency-independent) and 
factorization [lj| 

R(n) = i?23(ai,/3l, -OL\)R\2{ot2,p2, -"2) 

x R2 3 (a 3 ,0 3 ,-a 3 )e- i ^e- i ^, (11) 
hi =2C\\ — C22 — C33, fi2 = \{&22 — C33), Cij := (ij^aj, 
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with SU(2) subgroup matrices 
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R2s(a, (3, -a) = I cos § -e~ la sin f 
,0 c M sin§ cosf 



R 12 (a,/3,-a) = \ e ia sh 




immanant and determinant of and the _D-functions 

of irreps of SU(3). 

Using Young tableaux to denote the corresponding 
(12) functions of the matrix R(Q) constructed with frequen- 
cies of the output modes fixed with respect to those of 
the input modes, we observe the following. 







■ ( 13 ) Observations 1. For any R(Q) eSU(3) 
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We employ Young tableaux methods to determine out- 
put states given the three-photon input state 



I 111) = / dwi j duj 3 j dcj3(/>(wi)0(w 2 )<Kw 3 ) 

xe - 2 r le ^ 3 r 2 | 1(wi)1(w2)1(w3)) S (M) 

with time delays t\ and ti between modes 1 and 2 and 
between modes 1 and 3, respectively. The three-photon 
input state (fT4|) is 



□ 



□ 
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cm © EP © EP © § 



(1,0)® (1,0) 8(1,0) -> (3, 0)9(1,1) 9 (1,1) 9(0,0) 

(15) 

with the second row giving the standard labels in terms 
of highest weight. The last element of ([15")) is the fully 
antisymmetric irrep (0,0) with SU(3) labels (A, /it) such 
that A = 1 — 1=0 is the difference between the number 
of Lis in the first and second rows, and /U = 1 — 1 = 
is the difference between the number of Ds in the second 
and third rows. 

Basis states for the irrep (A, /i) are denoted 
|(A, \x)v\v?.vz\ I) with Vi the number of photons in port i 
and ^1 + ^2 + ^3 = A + 2/i. The index / distinguishes 
states with the same weight [y\ — V2, V2 — ^3) belonging 
to different irreps of the SU 2 3(2) subgroup of SU(3). In 
this notation, the three-photon state decomposes into 

|1( W1 )1M1M) = ^ |(00)000; 0) + ^ 1(30)111; 1) 
+ 1 1(11)111; 0) 1 + ^= |(11)111; IK 



I (11)111; 0) 2 



2- l(ii)iii; 1)2 



(16) 



The SU(3) irrep (1, 1) occurs twice in Eq. (fT5)l ; thus, an 
additional subscript (1 and 2) is needed to denote basis 
states that belong to these distinct copies of (1, 1). This 
input transforms as i?(fi)| I(a;i)l(aj2)l(w3)) and can be 
expanded in terms of the appropriate SU(3) /^-functions. 

The Young tableaux also label the representations cm, 
EP and S of 53 (the six-element permutation group of three 
objects). Characters of cm, EP and I are needed to con- 
struct the permanent, immanant and determinant of a 
3x3 matrix [13], respectively. Whereas the connection 
between coincidences at the interferometer output, the 
D- functions, and the permanents is clear [il[lll[lH, it be- 
hooves us to find a relationship between the permanent, 



Per(R(n)) EEcm(tt) = D^ )mn)1 (fl) 



Imm(R(Q)) =ff(fi) = , in)1 (fl) + D^'M 



(1,1) 



Det{R{tt)) eeB(Q) = D 



(in)l;(iil)i 

(0,0) 



(111)0;(1 



(000)0:(000) 



so the full permanent of the general 3x3 SU(3) matrix 
equals the D -function for irrep (3,0) with weight input 
and output. The immanant of R(Cl) is the sum of D- 
functions for the irrep (1,1) with weight input and out- 
put, and the determinant of R(fl) is just the D-function 
for irrep (0,0). 

Observation 2. For R(SY)kj the 2x2 submatrix of R(il) 
with the k th row and j column removed, 



PeriR(n) kj ) =nqy(fi) =d 



for ^(lioj.Uj.-mo- u.(^) the D function with in output 
state k and 1 otherwise, and in the input state j and 1 
otherwise. For instance 



(2,0) 

(liOfeKfcKiiOj 



Per(R(fl) 13 ) =mi 3 (0) 



D {2 ' 0) (CI) 

ly (011)l;(110)i^ 2 ' 1 - 



(17) 



In view of Observations 1, we consider measuring coin- 
cidences for the monochromatic (continuous- wave) three- 
photon input state (l(wi)l(w 2 )l(a;3)| s : 

S (l( Wl )l(^ 2 )l(w3)|i?(^)|lK)l(w2)l(^3)) S 



6- ty (lll)l;(lll)l\"'' 
-,(1,1) 



-D 



3 (^(lll)l;(lll)l(^) + ^(111)0;(111)0 



( °' 0) (fl) 
(000)0;(000)0V"-' 

(«)) 



= Wfi) + §EP(0) + il(fi). 



(18) 



As Eq. (|18|) is covariant under permutation of the output 
frequencies, detecting photons of frequencies (wi, u)j, Uk) 
in output ports (1, 2,3) respectively implies that 



D 



(A„u) 

(iii)J,(in) 



lllu (n) = <(A,M)(in)J|ii(n) 1(^^(111)/) 



h> <(A, /i)(lll)J| ( Pyfc )" 1 ^)l(A,M)(lll)I) (19) 



with pijk{coi, u>2, UJ3) H> (uii, uij, ujk) a permutation op- 
erator, and the action of pijk on SU(3) basis states is 
known (lij . 

Permuting output frequencies effects permutations of 
the rows in the basic SU(3) matrix. The permanent thus 
remains unchanged; the determinant picks up a sign for 
odd permutations; and the immanant transforms in a 
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more complicated way according to the 2-dimcnsional 
irrep of S3. Equation (|18l) remains valid provided the 
above changes are implemented. 

The coincidence rate for output modes (1,2,3) is 



Pm(f2) = s (m|i? t (ft)ni ® ni (g) nii?(fi)|m) 



which is a function of the two inter-photon delay times t\ 
and T2 and plotted in Fig. [I] for various values of ft and 
choices of the three detector spectral terms u)i and crj, 
i = 1,2, and 3. This rate Pm(f2) is a sum of six terms, 
one for each possible permutation of the three frequencies 
(uii, ujj 7 ujk) of the photons in output channels (1,2,3) re- 
spectively Each term is the product of an fi-dependent 
interferometer term as in Eq. (fT5|) and a detector spectral 
term as in Eq. $2$. 

The permanent is invariant under permutation of out- 
put frequencies: the coefficient of Z)( 3 '°) is constant for 
all output frequencies; upon rearrangement, Pm(fi) con- 
tains a single £K 3 >°) function multiplying the sum of all 
detector terms. Interferometer settings that make £)( 3 >°) 
vanish also make the permanent vanish from the coin- 
cidence rate, thereby resulting in a dip with respect to 
delay parameters Ti and T2 when detector terms are cho- 
sen to eliminate, by symmetry, the contributions from 
the determinant and the immanants. 

In Figs.ffia) and (b), we choose such that Z)( 3 > ) = 
and £) (3 ' ' = — l/4v% respectvely. Detector terms are 
partially symmetric under exchange 1 2, thus extin- 
guishing contributions from the immanants and determi- 
nant. As expected, only Fig. [TJa) has vanishing dip at 
the origin of the plot corresponding to two zero time de- 
lays. We also provide two examples of three-photon coin- 
cidence rates to demonstrate the complexity of dips that 
can arise when imperfect sources and detectors are used. 
In Figs. [Ifc) and (d), L>( 3 <°) = 0, but detector terms 
are asymmetric: contributions from the immanants and 
determinant interfere, and we obtain dips at non-zero 
time delays. Although here we have focused on the rate 
of obtaining three-fold coincidences corresponding to de- 
tecting one photon at each output port, we can easily 
calculate other coincidence rates such as measuring i>i 
photons exiting the first port, 1/2 from the second, and 1/3 
from the third port simply by replacing III ® III IT in 
the calculation of the rate Pm(f2) by Ii Vl <8 n„ 2 CED n„ 3 . 

In summary we have used SU(3) group theory to cal- 
culate the photon-coincidence rates at the output of a 
three-channel interferometer given single photons enter- 
ing each of the three input ports. We employ Young 
tableaux to decompose the direct product of irreps corre- 
sponding to each input photon into a direct sum of sym- 
metric, antisymmetric, and partially symmetric irreps. 
Our technique provides a powerful calculational tool for 
modeling and interpreting output data from realistic ex- 
periments that arc on the cusp of taking the two-photon 
two-channel Hong-Ou-Mandcl dip to a new and excit- 
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FIG. 1: A plot of the coincidence rate with for three dif- 
ferent sets of Q. = (qi, / 9i,Q2,/3 2 ,Q3,/93,7i,72): (a) fi = 
(0,7r/2,7r/2,7r,0,7r/2, tt/2, tt), uj = 0, a = 0.1, Wi = 0, 
<Ji = 0.1, 0J2 = 0, <T2 = 0.1, and W3 = 0, CT3 = 1, 

(b) n = (0,7172,0,7172,0,7172,0,0), ujq = 0, 0-0 = 1, wi = 0, 
<ti = 0.1, L02 = 0, <T2 = 0.1 and 0J3 = 0, as — 0.01, 

(c) fi = (0,7r/2,7r/2,2cos" 1 (l/v / 3),0,7r/2,0,0), too = 0, 
(jo = 0.5, uji = 3, ai = 0.2, 0J2 = 2, 02 = 0.2 and ujz = 1, 
(73 = 0.2, and (d) same Q. configuration as (a) but with 

0J = 0.1, (7 = 0.1, Wl = 0.95, (71 = 0.11, U) 2 = 0, (72 = 0.1 

and <J3 — 0, LJ3 — 0.99. 



ing regime of single photons entering multi-channel in- 
terferometers. Although brute- force calculations can be 
used to model the outputs from such interferometers, our 
methods are much more powerful than typical theoretical 
quantum optics techniques for modeling photonic quan- 
tum interferometry, both in terms of making the model 
tractable and also in making the resultant dips under- 
standable. We have provided a novel and clear connec- 
tion between the Wigner D functions for SU(3) and per- 
manents, immanants, and determinants, and this connec- 
tion is central to understanding the profiles of the photon 
coincidence dips as depicted in Fig. [T] This connection 
is especially important in light of the potential applica- 
tion of higher-order photon-coincidence dips to quantum 
computation such as for the boson sampling problem [Io| . 
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